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Abstract: The energy dependent Green's function {E — H)~^ for an interface Hamil- 
tonian which interpolates between two and three dimensions can be calculated explic- 
itly. This yields an expression for the density of states g{E, zq) on the interface which 
interpolates continuously between the two-dimensional g{E) =constant behavior for 
high energies and the three-dimensional q{E) oc ^/E behavior for low energies. 

PACS: 05.30.Fk (Fermion systems and electron gas), 71.10.Pm (Electrons in reduced 
dimensions), 73.20.-r (Electron states at surfaces and interfaces). 



1 Introduction 



Many quantities of physical interest depend on the number d of spatial dimensions. 
This includes potentials and two-point correlation functions, which are proportional 
to the zero energy Green's function 



G(.)(r) 



-r/2, 
-(27r)-Mn(r/a), 



d=l, 
d = 2, 

d>3, 



(1) 



as well as densities of states, which e.g. for non-relativistic free particles are 



g^,^(E)=ge{E)J- 




E 



■d-2 
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These are densities of states per c?- dimensional volume and per unit of energy. The 
factor g is the number of helicity or spin states of the particles. 
The corresponding dependence of the relation between the Fermi energy and the 
density n of electrons on d is 



n 



id) 



¥T{{d + 2)/2)\ 271 
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The generic case of physical interest is = 3, of course. Yet we frequently use mathe- 
matical techniques in d = 2 for the theoretical modeling of electrons or quasi-particles 
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on surfaces or interfaces. But how two-dimensional is e.g. an ensemble of electrons in 
an interface? Is the two-dimensional logarithmic two-point correlation function ap- 
propriate for the description of a gas of particles in an interface? Or should we rather 
expect the three-dimensional correlation? The correct answer will certainly lie 
somewhere in between, and how much in between will depend both on the specific 
system and its parameters. 

To address these kinds of questions analytically, dimensionally hybrid Hamiltonians 
of the form 



H = ^jd'^xjdz ( ■ Vip + d,^^ ■ d,ip) + Jd^^J dz ip^Vilj (4) 



2JI 



z=0 



were introduced in [T]. Here the convention is to use vector notation x = {x,y), 
V = {dx,dy) for directions parallel to the interface, while z is orthogonal to the in- 
terface. 

The potential term will generically also include three-dimensional and two-dimensional 
terms, 

V^(x,^) = H(x, z) + V,{^)6{z), 

but the competition between two-dimensional and three-dimensional behavior of 
physical quantities turns out to be mostly a consequence of competition between 
the two-dimensional and three-dimensional kinetic terms. 

The two-dimensional mass parameter is a mass per length. In simple models it is 
given by 

m 

where L± is a bulk penetration depth of states bound to the interface at 2; = 0, see 
Section El 

The zero energy Green's function for the Hamiltonian 



Ho = —Jd^x J dz (V^+ ■ + ■ d,ip) + ^ y rf^x ■ 



(5) 



2=0 



for perturbations in the interface {z' = 0, G(^x — x', z) = (x, ^IGlx', 0)) satisfies 



m 



(A + dl)G{^ - x', z) + —6{z)AG{x - x', 0) = -<5(x - ^')6{z) 

/i 

and was found in [Tj (r = |x — x'|), 

G{^-^,z) = -J^ dk^^^Mkr), i=-. (6) 
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Figure 1: The upper dotted line is the three-dimensional Green's function (47rr)~^ in 
units of the continuous line is the Green's function (^^, and the lower dotted line 
is the two-dimensional logarithmic Green's function, x = r/i. 



The Green's function in the interface is given in terms of a Struve function and a 
Bessel function, 



1 



G(x-x',0) = — 



and interpolates between two-dimensional and three-dimensional distance laws, 



r<£: G(x-x,0) 
r>£: G(x-x',0) 



1 

Anr 



-7 - ln(^ 



r 

2 / H 



(7) 

(8) 
(9) 



For the model ([5]) the Green's function G{'x, 0) appears in the first place as a two-point 
correlation function in the interface, but it can also be realized as an electromagnetic 
potential for an electromagnetic Hamiltonian 



H[F] = ^Jd'^Jdz + El + B^ + Bl) + £j rf^x (E^ + Bl) 



(10) 

2=0 



if the fields which are continuous across the interface yield a special interface contri- 
bution due to a finite limit 

2£ = lim CrL. 

L=0 

Here is the relative permittivity of the interface and L is its thickness. 
The correspondence between electronic two-point correlation functions and electro- 
static potentials is often used to map the quantum partition function of a free 
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fermion gas into the partition function of a classical Coulomb gas or plasma. The 
two-dimensionial Coulomb gas also plays an important role in the variational treat- 
ment of the fractional quantum Hall effect [21 [3l H]. The appearance of G(x, 0) as 
the Green's function for both of the dimensionally hybrid Hamiltonians ([5]) and ( fTOl) 
indicates that the free fermion gas - Coulomb gas duality persists in the transition 
regime between two and three dimensions. However, the objective of the present pa- 
per is the extension of G to non-zero energy and for perturbations off the interface, 
and the discussion of implications on the density of states. 

I will further elaborate on the motivation and justification for dimensionally hybrid 
Hamiltonians of the form (j4]) in Section O The calculation of the energy dependent 
Green's functions and the resulting density of states on the interface will be discussed 
in Sections [3] and HI Section [5] contains the conclusions. 



2 Dimensionally hybrid Hamiltonians 

The interesting part about the Hamiltonian (jlj) is the competition between two- 
dimensional and three-dimensional kinetic terms. Hamiltonians with low-dimensional 
6 contributions to the potential V are standard in quantum mechanics, and we easily 
understand the model assumptions behind the use of low-dimensional potentials. On 
the other hand, Hamiltonians with competing low- dimensional and three-dimensional 
kinetic terms are certainly not as common as superpositions of low-dimensional and 
three-dimensional potentials. I will denote Hamiltonians with superpositions of low- 
dimensional and three-dimensional kinetic terms as dimensionally hybrid Hamiltoni- 
ans. The introduction of these Hamiltonians in [1] was based on physical intuition 
and mathematical curiosity. The primary physical justification for the mixed kinetic 
terms is electrons or quasi-particles which can propagate with an effective mass m^, in 
a thin layer of thickness L and with mass m outside of the layer. In that case we would 
expect that particle propagation in the bulk-layer-bulk system should be described 
by a Hamiltonian of the form (jlj) with fi ~ m^^/L. The purpose of this section is to 
point out that models of the kind (jl]) can also be motivated in a different way from 
simple standard Hamiltonians with a standard three-dimensional kinetic term. 
As the simplest possible model consider an interface with an attractive potential 

l^(x, z) = -Wd{z), 

( 



The eigenvalues and eigenfunctions factorize with one set of states bound to the 
interface. 



/ ^ ^ / II ,\ exp(ik-x) VmW / m . 
^k,K,+ (x,2) = (x,z|k,ft,+) = — —exp\^-—W\z\ ) , (11, 
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2m ^ ^ 



2m 



m 
2^2 



W\ K 



mW 



and two orthogonal sets of transversally free states with transversal wavenumber 
k± > and energy E = h'^{'k'^ + k'j_)/2m, 



exp(ik ■ x) h'^k± cos{k±z) — mW sm{k_[_\z\) 



2tx 



/ ^ ^ / II 7 \ exp(ik-x) 1 . ^ 



27r -^TT 

The completeness relation for the factor in '?/'(x, 2;) = 0(x)x(2;) is [5] 

/■oo 
JO 

This implies a decomposition of the kinetic operator 

— y rf^xy rf^ (x+(^)V0+(x) ■ V0(x)x(^) + 5,x+(^)0+(x) ■ 0(x)9,x(;^)) 

- !L 
2m 



(12) 
(13) 



rf^x / dz (VV^+(x, ^) ■ VV^(x, z) + 9,^/^+(x, z) ■ dM^, z)) 



free 



+ |_ /"d2xV</.+ (x)-V</.(x) 
2m y 



2^2 



d X( 



where the first kinetic term on the right hand side only acts on components which 
are unbound in z direction, while the last two terms arise from the bound states. If 
we now express both the operators acting on free states and the operators acting on 
bound states in terms of ipi'x, z), we find 

fc2 r rco 

^Jd'^J dz (V^/'+(x, z) ■ VV^(x, z) + a,V^+(x, z) ■ dM^, ^)) 

fc2 r roo 

= ^Jd^^^J dz (V^+(x, z) ■ V^(x, z) + 9,^+(x, z) ■ 9,^(x, z)) 



free 



2mm 



rf2xV^/'+(x,0) ■ V^/'(x,0) 



bound 



w 
T 



d2x^/>+(x,0) ■ ^/'(x,0). 



This includes a superposition of kinetic terms similar to ([5]) with 

w?W m 
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where L±_ = is the bulk penetration depth of the bound states |k, n, +). 
Superposition of two-dimensional and three-dimensional kinetic terms might appear 
like an exotic concept for the description of particle propagation in the presence of 
surfaces or interfaces, but on a qualitative level the concept can be justified. And we 
have already seen in equations ([6]) and (JTj) that it allows us to provide simple analytic 
estimates on the impact of dimensional competition on propagation effects. In the 
following we will focus on a study of the competition between the different kinetic 
terms, i.e. we will neglect any potential terms and study the Hamiltonian Hq. 



3 The energy dependent Green's function 

The energy dependent Green's function Q{E) satisfies 
{E-H)g{E) = l. 

The Green's function G{E) which has ([6]) as the E = 0, z' = limit is related to g{E) 
through 

GiE) = 

Our objective is to find the energy dependent Green's function G{E) for the dimen- 
sionally hybrid Hamiltonian Hq in equation ([5]). For later comparison, we will first 
revisit the free Green's function Go{E) of the Hamiltonian Hq without the interface 
term, — > 0, and then calculate G{E). 

3.1 The free Green's function 

The free retarded Green's function is translation invariant, 
(x,z|Go(^)|x',z') = Go{E;^-x',z - z') 
and is well known to be 

Go(g;x,z) = ^^j^^^l^ e^p(-^s/-2mE{r^ + z^)\ (14) 

where we continue to use cylinder coordinates. For positive energy this corresponds 
to the standard choice of poles in Gq{E; k, k±) to generate outgoing spherical waves 
without any incoming spherical component, in symbolic notation 

^^^^ = E-H + ie =^ E-E^ + ie ^''^ 
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where u is a. degeneracy index. 

Green's functions for surfaces or interfaces are commonly parametrized in an axially 
symmetric mixed representation like G{E] k, z, z'), see e.g. [6]. In bra-ket notation this 
corresponds for the free Green's function Go{E), which is also translation invariant 
in z direction, to 

{k, z\Go{E)\k' , z') = Go{E;k,z-z')6{k-k'). 

For later comparison we will also briefly recall the explicit form of the free Green's 
function Gq{E) in the axially symmetric mixed parametrization. The equation 

(5.'-k^ + ^) G,{E,k,z) = -5{z) 
yields 

^ / X 1 /" expfi/c 1 2) 

Gq{E: k,z) = — / dk^ , , ^, J,^,, 16 

' ' ^ 2n J kl + k^- {2mE/h?) - ie ^ ' 

neih^k^ - 2mE) ( r—— -\z\\ 

= : ^ pyp -Vh^k^ - 2mE^-^ 

2Vh?k^ - 2mE V ^ J 

2V2mE - W-k^ V ^ / 

Note that this axially symmetric representation of the free Green's function depends 
on the number d — \ oi dimensions perpendicular to z only through the number of 
components of the wavevector k. 



3.2 The interface Green's function G{E) 

Our objective is to generalize the Green's function G(x — x', = (x, z|G'(0)|x', 0) ^ 
both to general values of E and z' , (x, 2;|G(ii^)|x', z'). We will also keep the general 
value ^0 for the location of the interface and investigate the dependence of GiE) on 
this parameter. 
The Hamiltonian 

iJo = y rf'x y dz (V^A+ ■ VV^ + ■ d^i)) + - ^o) VV^+ ■ VV^^ 

yields the Schrodinger equation 

E^A(x, z) = -^{A + d!)^lj{^, z) - ^6{z - zo)A^A(x, z), 



or 
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The last equation in (x, z) representation is 
Substitution of the Fourier transform 



—E + A + 9f + S{z - zo)-A (x, z\GiE)\^', z') = -5(x - x')5(^ - z'). (17) 



(x,z|G(E)|x',z') = ^ jd^^ jd^k' {k,z\G{E)\k',z')exp[i{k-^-k' ■^')] 
yields 

0,771 \ 771 

-^E~k' + dlj (k, z\G{E)\k\ z') - -k'6{z - zo){k, z\G{E)\k\ z') (18) 
= -6{k-k')6{z- z'). 
This yields with 

{k, z\G{E)\k' , z') = {z\G{E,k)\z')6(k-k') 
the condition 

(2777 \ 777 

—E -k^ + dlj {z\G{E, k)\z') - -kH{z - z^){z\G{E, k)\z') = -5{z - z'). 

Fourier transformation with respect to z yields 

'2777 \ 

—E~k^-kl\ {k^\G{E,k)\z') (19) 



-k^ / (iK_L exp[i(K_L — k^)zQ\{K^\G{E,k)\z') = ^= exp(— ifc_L2;') 

' J V27r 



This result implies that {k±\G{E,k)\z') has the form 

(exp[iA;^(zo - ^')]/v^) + f{E, k, z') 



eM^k,z,){k,\GiE,k)W) ^ _ (2^^/^,) 

with the yet to be determined function f{E, k, z') satisfying 

For the treatment of the integrals we should be consistent with the calculation of the 
free retarded Green's function (fT6!) . 



/ ^'^^ exp(i^±^) _ ^o/^2i 2 



i^e{2mE - h^k'' 



exp(—y/h?k'^ — 


2mE 


\z\/h) 


Vh?k^ - 

exp(iy/2mE — 


2mE 

h^k^l 


z\/h) 


V2mE - 


- h^k^ 





This yields 

H k 

mh 



eifi^k^ - 2mE) , ^e{2mE - h^k^' 



2/iV27r 



expl — V h?'k'^ — 2mE 



V¥W^^2rnE V2mE - h?k^ 
e{h'^k'^ - 2mE) 

Vh^k^ - 2^ 

^/2mE - h?k^ V ^ 



fiE,k,z') 



h 



and therefore 
{k^\G{E,k)\z) 



1 



y^kl + k^ -{2mE/n?) 
hkHe{!fk'^ - 2mE) 



exp(— i/c_Lz') 



(20) 



^h?k^ - 2mE + nk2£ 

. ;ik2£e(2mE-;i2i^2) 

~ \2mE -rfi^ + \hk^E 
Here we use the definitioE0 

£= — 
2/i' 

Fourier transformation with respect to k_[_ yields finally 

nein^k^ - 2mE) 



expl — iA;_L-2o — VK^k"^ — 2mE- 



exp ( —ik±zo + W2mE — ffik^- 



{z\G{E,k)\z') 



2Vh?k'^ - 2mE 
hkH 



expl — V/l^k^ — 2mE 



\z — z 

h 



(21) 



^Jh^k^ - 2mE + hk^^ 



expl — V fi^k^ — 2m£' 



z — Zo I + I ^' — ^0 I 



.ne(2mE-n2k2^ 



+ 1 



2^2mE - ;i2k2 
;ik2£ 



z — z 



V2mE - h^k^ + i;ik2£ V 



exp ( iV2mE — h^k^ 
V ^ 

exp ( W2mE — h'^k'^- 



h 



^ol + F 



2^01 



h 



This result is translation invariant in the transverse z direction for scattering off 
perturbations on the interface, z' = zq, 



{z\GiE,k)\z') 



G{E,k,z-zo). 



^ In the simple model from Section[2l £ ~ L±/2 = (2k)~^ would be the bulk penetration depth 
of the probability densities | (x, z|k, k, +) p of the bound states. But note that we have neglected any 
potential contribution in Hq, such that the results derived here are not directly applicable to the 
model from Section O 
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4 Density of states and Fermi energy on the inter- 
face in the dimensionally hybrid model 



The equation ( 1151) yields a standard expression for the density of states in terms of 
the imaginary part of Green's functions, see e.g. [7], 



Q{En,y^,z) = 5( V'(x,2;|n, z/)(n, i/|x,2) = S(x, z|^(£'„)|x, 2;) (22) 

z — / 7J- 

V 

2mg^ 



■53(x,z|G'(E„)|x, z). 



Here we explicitly included a factor g for the number of spin or helicity states, because 
the summation over degeneracy indices in f[T^ usually only involves orbital indices. 
For translation invariant Green's functions 

{:>c,z\Q{E)\x',z') = ^(i?;x-x,2;- z') 

we have 

g{E) = ~^^g{E- X = 0, z = 0) = --^^1^55 j rf-^^^k k, z = 0). (23) 

Insertion of the free retarded propagator (fT6l) reproduces the standard density of 
states ([2]), of course. 



2mg he(E) rV2^/n 



Trh^ 2^-iv^''"'r(((i- l)/2) Jo V2mE - h^k^ 

^ , /m \fE 
gQ{E). ' 




2-K T{d/2)¥' 



The interface at Zq breaks translational invariance in z direction, and we have with 
{k,z\G{E)\k',z) = {z\G{E,k)\z)6{k~k') 

q{E,z) = ^53(x,;.|G(E)|x,.) = J {z\G{E,k)\z)- 

We will use the result ( 12T]) to calculate the density of states g{E, zq) for the Hamil- 
tonian ([5]) on the interface. Substitution yields 



0{E,zo) = j£k{zo\GiE,k\zo) 



^^'^ „ , ^2mE - h^k^ 



gm 

2^^'"'7n ^'"'"2mE -^i^k^ + h^kH^' 



e{E) / dkk 
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Figure 2: The dotted line is the three-dimensional density of states in units of 
gm/Air'^h'^i. The continuous line is the density of states fl2^ . < x = 8mE£'^/h? < 1. 



The evaluation of the integral yields 
gme{E) 



Anmi^n? - 8mEP 



Q{h' - 8mEt 



(24) 



k+VW^S^) arctanf ] 
- (h-Vh-'-SmEiA arctanf ^^^^^ 



gmQ{8mEe^ - 



h 



VSmEi'^ - ^2 



In 



+ arctan 



V8mE& + iVSrnE 



iVSrnE - VSmEP - It? 

h 

arctan 



This is a more complicated result than the density ([2]) for d = 2 or d = 3. However, 
it reduces to either the two-dimensional or three-dimensional density of states in the 
appropriate limits. For large energies, i.e. if the states only probe length scales smaller 
than the transverse penetration depth £, we find the two-dimensional density of states 
properly rescaled by a dimensional factor to reflect that it is a density of states per 
three-dimensional volume. 



8mEf > : g{E, Zq) e{E) 



gm 



-r-.Q{d=2){E). 



(25) 
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For small energies, i.e. if the states probe length scales larger than the transverse 
penetration depth we find the three-dimensional density of states 



gym 



E = g^a=3)iE). 



(26) 



This limiting behavior for interpolation between two and three dimensions is con- 
sistent with what was observed already for the zero energy Green's function in the 
interface ([TIE]). 

Equation (l2ll) also implies interpolating behavior for the relation between electron 
density and Fermi energy on the interface for the Hamiltonian ([5]). The full relation 
is with g = 2 



n{zo) 



Ep 



+ 



+ 
+ 
+ 



dE g{E, zq) 
^/mEp 1 

0(^2 - SmEpf) 

AmEpf + 
QiSmEpf - h^) 



(27) 



AmEpf - h\/ h? - SmEpP ) arctan 



8mEp 



arctan 



AmE-^^ 



871^ hp 

F 



^SmEpP - In 



^/SmEpl - VSmEpP - 

h 



h 

AmEp"^ 



arctan 



fVSmEpi + VSmEpP - h? 

h 



-arctan 



/VSmEpi - V8mEpP - 

' h 



This approximates two-dimensional behavior for niEpP ^ h"^, 
, , mEp 1 



-n 



(d=2) , 



and three-dimensional behavior for mEpP' -C /i^, 

-3 

= '^(d=3)- 



77,(2:0) 



^2mEF 



It is intuitively understandable that the presence of a layer increases the Fermi energy 
for a given density of electrons. The presence of a layer implies boundary or matching 
conditions which reduce the number of available states at a given energy. 
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5 Conclusion and Outlook 



Inclusion of competing kinetic terms for effective propagation of particles in low- 
dimensional subsystems and bulk materials implies dimensional interpolation effects 
for Green's functions and for quantities derived from the Green's functions in such 
systems. Simple model systems with competing kinetic terms can be solved exactly 
and yield analytic insights into the transition between two-dimensional and three- 
dimensional properties in materials with low-dimensional subsystems. 
The present paper focused on the study of the competition between two-dimensional 
and three-dimensional kinetic terms, thereby neglecting any bulk and layer potentials. 
In that case the transition scale between two-dimensional and three-dimensional be- 
havior is given by £ ~ m/yU ~ {m/m^)L. Please note that inclusion of the confining 
layer potentials will change this scale for most systems. The approximation Hq per se 
as an approximation for specific systems is only useful for the study of propagation 
effects and impurity scattering of particles which are not strongly bound to a thin 
layer, but which are affected by its presence to the extent that propagation in the 
thin layer is described by an effective mass. 

The purpose of the present investigation was to further advance a novel tool for the 
study of low- dimensional systems, not to derive generic quantitative properties of 
these systems. The Green's function G{E) calculated here should nevertheless prove 
useful for the study of impurity scattering of weakly coupled particles in thin layer 
systems, with the impurity potentials treated as perturbations. 

Acknowledgement: This work was supported in part by NSERC Canada. 

References 



[1] R. Dick, Int. J. Theor. Phys. 42 (2003) 569-581, cond-mat/0204534 



[2] R.B. Laughlin, Phys. Rev. Lett. 50 (1983) 1395-1398. 

[3] T. Chakraborty, P. Pietilainen, The Quantum Hall Effects, 2nd ed., Springer- 
Verlag, Berlin 1995. 

[4] T. Chakraborty, Adv. Phys. 49 (2000) 959-1014. 

[5] S.H. Patil, Am. J. Phys. 68 (2000) 712-714. 

[6] M. Lannoo, P. Friedel, Atomic and Electronic Structure of Surfaces, Springer- 
Verlag, Berlin 1991. 

[7] S. Doniach, E.H. Sondheimer, Green's Functions for Solid State Physicists, Im- 
perial College Press, London 1998. 



13 



